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The effects of an electron-phonon (e-ph) interaction on the thermoelectric properties of Na^CoC^ 
are analyzed. By means of dynamical mean field theory calculations we find that the e-ph coupling 
acts in a cooperative way with the disorder, enhancing the effective binary disorder potential strength 
on the Co sites, which stems from the presence or absence of a neighboring Na atom. Hence, the 
inclusion of the e-ph coupling allows us to assume smaller values of the binary disorder potential 
strength - which for Na^OoC^ are in fact also more reasonable. More generally, we can conclude 
that the interplay between disorder effects and coupling to the lattice can be exploited to engineer 
more efficient thermoelectric materials. 



Materials with a high thermopower— attract strong in- 
terest because of their potential use in heat-electricity 
conversion. Recently, experimental and theoretical in- 
vestigations evidenced that materials with strongly cor- 
related electrons 2 - can have very high values of the ther- 
moelectric figures of merit ZT, i.e., a high thermoelectric 
efficiency In this respect, Na x Co02 is one of the most 
promising materials^ besides presently used semiconduc- 
tors such as PbTe4 It has been shown that several condi- 
tions are important for getting the high thermopower: a 
"pudding mold" -like shape of the band structure^, strong 
correlation and also disorder— ~—. 

Phonons have also been argued to have a positive ef- 
fect on the thermopowe r 10 ' 11 : "nano phonon engineer- 
ing" as in heterostructures, thin films, cage structures 
and nanowires can indeed overcome the obstacle that 
a high electrical conductivity is usually connected with 
a high thermal conductivity, which is inversely propor- 
tional to the figure of merit. However, in materials char- 
acterized by the presence of several concurrent physical 
mechanisms, not only the individual effects, but also the 
interplay and the possible synergies between them have 
to be carefully analyzed. For instance the cooperation 
between the electron-phonon interaction and the disor- 
der is in this respect highly unexplored^. 

The importance of the Na disorder potential was first 
brought up by Singh and Kasinathan^ who attributed 
the differences between their local density approxima- 
tion (LDA) and photoemission experiments^ to the dis- 
order. The effect of two inequivalent sites on the elec- 
tronic correlations was then analyzed by Marianetti and 
Kotliar^ii^ by means of LDA+dynamical mean field the- 
ory (DMFT)2i with the inclusion of the disorder at level 
of the coherent potential approximation (CPA). Further 
LDA+DMFT studies which focused on reproducing the 
thermopower of Na^CoO^^ analyzed the mechanism be- 
hind the influence of disorder on electronic correlations: 
because of the disorder the lattice sites have locally dif- 
ferent occupations and since a fraction of them moves 



towards half-filling, correlation effects get effectively en- 
hanced. This increase in correlation leads to stronger 
electron-hole asymmetry in the LDA+DMFT spectra 
w.r.t. the LDA and this favors a large thermopower. 

An important point of the LDA+DMFT studies dis- 
cussed above is that the value used for the disorder po- 
tential Ae = 0.55eV is about ten times larger than the 
one which can be estimated from the energy distance of 
the LDA bands attributed to the inequivalent sitesii It 
is therefore evident that electronic correlation and dis- 
order effects arc not enough to explain the experimental 
findings. This leaves the question of which are the rele- 
vant interactions to get high values for the figure of merit 
answered. 

Here we propose the cooperative interplay between 
electron-phonon interaction and disorder as a new mecha- 
nism for increasing the thermopower. Our findings there- 
fore suggest a very simple route to high figures of merit 
in strongly correlated materials with a favorable shape of 
the band structure: Introduce a small amount of disorder 
in those compounds in which the electron-phonon inter- 
action is also important. To support our claim we add 
an electron-phonon coupling to the LDA+DMFT calcu- 
lations for Na^OoC^, i.e. to the simplest thermoelectric 
materials with electronic correlation, disorder and favor- 
able shape of the bands around the Fermi level. We fo- 
cus on the characteristics of the self-energy that in Refs. 
0H were identified to enhance the thermopower. The 
goal is to obtain the same features in the presence of an 
electron-phonon coupling, without the need of assuming 
unrealistically large values of the disorder strength. We 
emphasize that the resulting phonon enhancement of the 
thermopower which is introduced here is an additional ef- 
fect besides the suppression of the thermal conductivity 
discussed befor o 10 i n . 

Before we present our results for the specific case of 
Na^OoC^, let us discuss why we expect in general e- 
ph interaction and disorder to pull together in the same 
direction: A local e-ph coupling to the electronic charge 
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gives - among other effects - a shift of the local energy 
level, i.e. an effective shift of the local chemical potential. 
The role of such a shift has been investigated by one of 
us before^ in a more general context of DMFT model 
calculations: For small values of the phonon frequency it 
is well described by the following linear behavior: 

u(A) ~ fj, - nXW (1) 

where n is the density, W the electronic bandwidth and A 
the dimensionless electron-phonon coupling (see below). 
This means that the e-ph coupling A reduces the local 
chemical potential /i(A) of the occupied sites (n > 0). In 
the CPA, i.e. the framework used here to treat binary 
disorder in DMFT, a local potential Ae is added to the 
sites where disorder is present. Sites with and without 
disorder have hence different occupations. Therefore the 
e-ph coupling will act quite differently on them and will 
"amplify" the difference in the level renormalization be- 
tween the two sites. This simple mechanism results in an 
enhancement of the disorder effects, as we will argue and 
also demonstrate numerically below. 

We employ a single ai g band model for Na^OoC^ de- 
rived from the LDA, supplemented by a local Coulomb 
interaction U = 3.5 eV: 

H e = ~ l ^ c % c ia + U Y1 n n n u + Ae ^ n; (2) 

ij.fr i iGVac 

Here, c^,. and Ci a create and annihilate an electron on site 
i with spin a, respectively, n ia c^ a c ia is the electron 
density on site i for the spin a and denotes the electron 
number operator on site i summed over spin. As already 
introduced in Ref. [H, sites with an adjacent Na atom 
(i G Na) arc not altered, while a local disorder potential 
Ae is added to the "vacuum" sites (i £ Vac). 

The coupling to the Ai g phonon is described by the 
following Holstein form: 

H e _ ph = -g 2J (a+ + en) ma + ljq ^ a+a- . (3) 

i.a i 

The clcctron-phonon coupling is given by A = Eu po i/W 
with Eupoi = 2g 2 /wo being the bipolaronic binding en- 
ergy. This Hamiltonian captures the local part of the 
electron-phonon interaction in Na K Co02, which as shown 
in Ref. |23j, is the one with the largest coupling. The ad- 
ditional physics associated to the non-local part of the 
electron-phonon coupling lies beyond the scope of the 
DMFT approach which we will adopt here. 

The impurity solver used for the DMFT equations is a 
continuous-time quantum Monte Carlo in the "CT-INT" 
implementation which can take into account the phonon 
degrees of freedom. This algorithm is very well suited to 
include a local Holstein electron-phonon coupling to an 
Einstein mode. The idea is to integrate out the bosonic 
mode so as to obtain a retarded density-density interac- 
tion. This retarded interaction can readily be taken into 
account within the framework of the CT-INT algorithm. 



Details of the implementation can be found in Ref. [24] 
and selected applications in Refs. I25II20I The CPA used 
for the disorder boils down to an arithmetic average of 
the local Green function weighted by the corresponding 
Na and vacuum densities riNa and nyac • For a given value 
of the disorder potential we adjust the chemical potential 
(i such that the total density n to t = a^Na + (1 — x)ny ac 
equals 1.7. 

In the absence of e-ph coupling Ae leads to a lo- 
cal charge density disproportionation An = ?iNa — ^Vac 
which, in the considered parameter range, turns out to 
be linear: Ae = 6 An. This is shown in Fig. [1] (red cir- 
cles). Using Eq. [T] and for a fixed charge imbalance, 
we can deduce within the framework of the CPA that 
Ae(A) = Ae(A = 0) - AnXW, thus leading to a A depen- 
dence of the coefficient b(X) = b — XW. Hence, in the 
presence of the e-ph interaction we conjecture that 

Ae = (6 - XW) An (4) 

In other words, for the same given disorder potential 
Ae, a larger charge disproportionation An is obtained if 
the e-ph coupling A is included. This explains the change 
in the slope upon increasing A observed in Fig. [TJ where 
the results of our calculation are shown. As we expect 
from Eq. (J4| the behavior is linear with and without e-ph 
coupling: As it should be, all linear fits (dashed lines in 
Fig. [JJ to the data have an intercept at (0, 0) and the lin- 
ear form accounts well for the numerical data. The com- 
parison with Eq. [4] is therefore qualitatively very good, 
albeit the actual values obtained for the slopes are larger 
than b — XW. Even though Eq. [JJ gives a quantitative 
estimate of the chemical potential, the compressibility, 
which corresponds to the inverse slope with respect to n, 
comes out less accurately. The interplay between elec- 
tronic and lattice interactions indeed may renormalize 
an effective electron phonon coupling, going beyond the 
approximation in Eq. [TJ 
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FIG. 1: Ae as a function of An for different values of A. 

The changes in the slope of Ae vs An visible in Fig. 
UJ can be also understood as follows: The compressibil- 
ity k = dn/dfi of a Fermi liquid in the presence of an 
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electron-phonon coupling read o 27 ' 28 : 

Ku (0) 



Ku(X) 



1 - ku(0)XW 



(5) 



where ku(0) is the compressibility in the absence of 
electron-phonon interaction. The Hubbard U enhances 
the compressibility w.r.t. the LDA value (in other words 
Ku(0) is bigger than the free compressibility given by the 
value of the density of the states) . The e-ph interaction 
then further increases it according to Eq. [SJ A qualita- 
tively similar increase can be observed in the slope of An 
as a function of Ae (actually a decrease of the inverse 
quantity, i.e. of the slope of Ae as a function of An), as 
shown in Fig. [TJ Since this is not the charge response 
to a uniform change of the chemical potential, we are 
not dealing with a true compressibility. An corresponds 
rather to a "CPA" -like shift of the relative occupations 
but the behavior with A is qualitative very similar to that 
of kjj (A) of Eq. [TJ Another difference is also that the true 
compressibility would not look as linear in fi as the CPA 
result for An as a function of Ae. Concerning the ex- 
act values for the slope with and without phonons, let us 
note that in previous model calculations within DMFT 
the compressibility was also only qualitatively described 
by Eq. [5] despite the fact that the chemical potential was 
very accurately given by the relation Eq. [if 2 ' 29 . 

We have therefore shown that the e-ph interaction ef- 
fects on the effective local potential go in the same direc- 
tion as those of the disorder. We have understood this in 
terms of the effect of an e-ph coupling on the compress- 
ibility of a Fermi liquid. For the purposes of the present 
analysis, this observation means that we can expect to 
obtain similar disorder-induced effects on the self-energy 
by considering the e-ph coupling and smaller (as well as 
more realistic) values of Ae compared to those used in 
calculations at A = 0. Thus, we can conclude that large 
values of the thermopower can be obtained in the pres- 
ence of e-ph interaction by combining a moderate degree 
of disorder with electronic correlation. 

To see this quantitatively we can concentrate on the 
signatures characterizing the large thermopower in the 
self-energy of Na^OoO^, namely - as shown in Refs. 0@ 
- a non-zero Im E(0) and a deviation from a linear- 
dependence of Im £(iu; n ) at small uj n , characterizing 
Fermi-liquid solutions. In Fig. [5J [3] and 2] we show the 
self-energy of the whole system ("Tot") and that of the 
individual sites with and without Na. Lines represent the 
case with no phonons and large strength of the disorder 
while the symbols indicates the cases with increasing A 
and correspondingly smaller and smaller Ae. From Fig. 
[2l [3] and [4] it is clear that in the presence of e-ph inter- 
action both self-energy features discussed above appear 
in the "Tot" curve for smaller values of Ae than without 
e-ph interaction. For A = 0.25, for example, a value of 
0.30 instead of 0.55 suffices (Fig. [2]) and the effect gets 
rapidly stronger upon further increasing A. 

It is also interesting to analyze how the individual 
self-energies are influenced by the electron phonon in- 




FIG. 2: Imaginary part of the self-energy on the Matsubara 
axis for electron-phonon coupling A = 0.25 (symbols) and A = 
0.00 (solid lines). 
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FIG. 3: Imaginary part of the self-energy on the Matsubara 
axis for A = 0.50 (symbols) and A = 0.00 (solid lines). 




FIG. 4: Imaginary part of the self-energy on the Matsubara 
axis for A = 0.75 (symbols) and A = 0.00 (solid lines). 
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teraction: The "Vac" self-energy, corresponding to a lo- 
cally higher density of holes is less strongly modified by 
phonons at small Matsubara frequencies, the reason for 
this being that the Hubbard interaction more effectively 
blocks charge fluctuations at higher densities. The "Na" 
self-energy does get instead more influenced by A, reflect- 
ing the fact the the "Na" sites are locally closer to the 
low-density hole limit in which the Hubbard U effects 
vanish. In this light, it is possible to also understand the 
different effectiveness of the e-ph coupling for Na^OoC^ 
and CPA disorder: the weight of the "Na" component 
is large and consequently, a modification of that compo- 
nent, as the one induced by A, results in a strong effect for 
the entire system. The picture of two components having 
locally two different site occupations explains also why 
the slope of the corresponding self-energies is changed in 
opposite ways: That of the "Na" self-energy increases in 
absolute value, while that of the "Vac" sites decreases. 
In the former the effects of U are indeed smaller and 
therefore the phonons act in a more conventional way, 
i.e. they tend to localize electrons (eventually leading to 
polaron trapping). For the "Vac" component, instead, 
the main effect of the e-ph interaction is a reduction of 
the total repulsion of the system. In other words phonons 
increase the distance from the metal-insulator transition 
U values and, as a result, electrons have a smaller ef- 
fective mass (smaller absolute value of the slope of Im 
E(iw n ))2fl. Both effects can be clearly observed in Fig. 
EUandffl 

Our analysis of the e-ph support to the thermopower 
of strongly correlated thermoelectric materials with dis- 
order, such as Na x Co02, would be a mere academic ex- 
ercise if the e-ph coupling would not be strong enough to 
give rise to a visible effect. In the literature there are two 
independent estimates of the e-ph coupling in this ma- 



terial. The first one is a direct LDA calculation within 
linear-response, where the values of A of 0.21 and 0.29 
are reported, for x = 0.3 and 0.7 respectively^. The sec- 
ond is an estimate based on electrostatic arguments and 
a comparison with Raman experiments^. This, for the 
(orbital) diagonal part of the A\ g oxygen mode gives a 
value of about 0.044, which summed over all q vectors 
gives about 0.13. A factor of two between the two num- 
bers is understandable because first of all the latter is 
the coupling to just one selected mode and also because 
it contains a very delicate estimate of the dielectric con- 
stant in Na a; Co02- In addition, the A estimated in LDA 
is extracted from the Eliashberg function and contains 
the phonon frequency renormalization effects, while the 
A estimated by Donkov, et al. is the "bare" e-ph coupling 
constant (for a discussion about the difference between 
these types of estimates see Ref. l32l ). 

If in our scheme we use a value close to the LDA es- 
timate for Na.j;Co02, namely A = 0.25, we get (Fig. [5]) 
strong effects on the self-energy for Ae = 0.30, i.e. almost 
a factor of two smaller than what used in calculations 
without e-ph interaction. Our analysis therefore fully 
supports the claim that there are the conditions to ex- 
ploit the e-ph interaction for enhancing the thermopower 
of strongly correlated materials, provided that the con- 
ditions of a favorable asymmetric shape of the bands is 
present and a moderate degree of disorder is included. 
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